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Abstract 

We present infrared spectroscopy of the Antennae Gal- 
axies (NGC 4038/4039) with NIRSPEC at the W. M. 
Keck Observatory. We imaged the star clusters in the 
vicinity of the southern nucleus (NGC 4039) in 0".39 see- 
ing in K-band using NIRSPEC 's slit- viewing camera. The 
brightest star cluster revealed in the near-IR (Mk(0) ~ 
— 17.9) is insignificant optically, but coincident with the 
highest surface brightne ss peak in the mid-IR (12 — 18/im) 
ISO image presented by Mirabel et al. (1998) . We obtained 
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Obscured star cluster 
and NIRSPEC slit 



NIESPEC SCAM 




high signal-to-noise 2.03—2.45 /j,m spectra of the nucleus 
and the obscured star cluster at R ~ 1900. 

The cluster is very young (age ~ 4 Myr), massive (M 
- 16 x 10 6 M Q ), and compact (density ~ 115 M pc~ 3 
within a 32 pc half-light radius) , assuming a Salpeter IMF 
(0.1—100 Mq). Its hot stars have a radiation field charac- 
terized by T G ff ~ 39, 000 K, and they ionize a compact H n 
region with n e ~ 10 4 cm~ 3 . The stars are deeply embed- 
ded in gas and dust (Ay ~ 9 — 10 mag), and their strong 
FUV field powers a clumpy photodissociation region with 
densities nn 10 5 cm~ 3 on scales of ~ 200 pc, radiating 
Lh 2 i-o s(i) = 9600 L . 

Key words: Galaxies: individual (NGC4038/39, Antennae 
Galaxies) - Galaxies: star clusters - Galaxies: H n regions 
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Figure 1. HST image of the Antennae overlaid with ISO 
12-18 urn contours ^Mirabel et al. 199$ ), and NIRSPEC 
K band image showing the slit position and several star 
clusters in the field. 



The Infrared Space Observatory (ISO) 12 — 18 (im im- 
age showed that the hot dust distribution is similar to 
that of the gas, but peaks at an otherwise incon spicuous 
point on th e s outhern edge of the overlap region ( Mirabel 



et al. 1998 ; see Wilson 200C ). This powerful starb urst knot 



1. Introduction 



is also a flat-spectrum radio continuum source ( Hummel 



The Antennae (NGC 4038/4039) are a pair of disk galaxies 
in an early stage of merging which contain numerous mas- 
sive super star clusters (SSCs) along their spiral arms and 
around their interaction region (Whitmore & Schweizer 



fc van der Hulst 1986) and may b e associated with an 
X-ray source ( Fabbiano et al. 1997 ). We imaged the re- 



1995; Whitmore et al. 1999). The molecular gas distribu- 



tion peaks at both nuclei and in the overlap region (Stan 
ford et al. 199C), but the gas is not yet undergoing a global 
starburst typical of more advanced mergers (Nikola et al 
1998). Star formation in starbursts appears to occur pref- 



erentially in SSCs. We chose to observe the Antennae be- 
cause their proximity permits an unusually detailed view 
of the first generation of merger-induced SSCs and their 
influence on the surrounding interstellar medium. 



gion around this knot, and discovered a bright compact 
star cluster coincident with the mid-IR peak. We obtained 
moderate-resolution (R ~ 1900) K-band spectra of both 
the obscured cluster and the NGC 4039 nucleus. 



2. Observations & Data Reduction 

NIRSPEC is a new facility infrared (0.95 — 5.6/im) spec- 
trometer for the Keck-II telescope, commissioned during 

1998|). It has a 



April through July, 1999 ( McLean et al. 
cross-dispersed cryogenic echelle with R ~ 25,000, and a 
low resolution mode with R ~ 2000. The spectrometer de- 
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Figure 2. NIRSPEC spectrum of the obscured star cluster 
shows nebular and fluorescent Hi emission with a contin- 
uum rising toward the red. Scaled sky counts are plotted at 
0.1 mjy. co-shaped curve represents an atmospheric COi 
band at 2.05 fim. 



Figure 3. NIRSPEC spectrum of NGC 4039 nucleus shows 
extended collisionally excited H2 emission and a strong 
stellar continuum marked by photospheric absorption. No 
Brj is present. Scaled sky counts are shown at 1.5 mJy. 



tector is a 1024 x 1024 InSb ALADDIN focal plane array, 
and the IR slit-viewing camera detector is a 256 x 256 
HgCdTe PICNIC array. 

We observed the Antennae with NIRSPEC during the 
June 1999 commissioning run. Figure I shows the HST 
and ISO image of [Mirabel et al. (1998)] together with NIR- 
SPEC slit-viewing camera (SCAM) images at 2 /zm; the 
SCAM images reveal that the mid-IR ISO peak is a bright 
(K = 14.6) compact star cluster located 20". 4 east and 
4". 7 north of the K-band nucleus. This cluster is associ- 
ated wiiJa » fain! (v = ^ Sj red (V-T = 9 Q) smuxe (# 80 



in Whjfcmace Xr Schsreizsr 1QQK) irigihlp with Spare TpIp. 



scope (Whitmore & Zhang, private communication). We 
obtained low resolution (R ~ 1900) A ~ 2.03 - 2.45 /jm 
spectra through a 0".57 x 42" slit at PA=77° located on 
the obscured star cluster and the nucleus of NGC 4039, 
for a total integration time of 2100 s. We reduced the data 
followi ng the standard proce dure, described in more de- 



tail by |Gilbert et al. (2000) , Reduced spectra are shown 



in Figures |2j and 



3. Massive Star Cluster 

The cluster spectrum is characterized by strong emission 
lines[] and a continuum (detected with SNR ~ 15) domi- 
nated by the light of hot, blue stars and dust. The nebular 
emission lines are slightly more extended than the contin- 



1 A table of measured line fluxes is available electronically 
from http: / /astro. berkeley.edu/~agilbert / antennae. 



uum, and the H2 emission is even more extended. This 
suggests a picture in which hot stars and dust are embed- 
ded in a giant compact H 11 region surrounded by clumpy 
(see §3.2) clouds of obscuring gas and dust whose surfaces 
are ionized and photodissociated by FUV photons escap- 
ing from the star cluster. 

For a distance to the Antennae of 19 Mpc (H p— 75 km 



s -1 Mpc" 1 , 1"= 93 pc) ( |Whitmorc et al. 1999| ), we find 



that the cluster has Mr = —16.8. We estimate the screen 
extinction to the cluster by assuming a ran ge of (V— K)p ^ 
— 1 as expected from Starburst 99 models ( Leitherer et al 
1999D , and that A K = 0.11 A v ( |Rieke fc Lebofsky 1985 1 ) . 
We find A v = 9 - 10 mag, which implies M K (0) = -17.9, 
adopting Ak = 1.1 (which is confirmed by our analysis 
of the H 11 recombination lines in §3.1). We can use the 
intrinsic brightness along with the Lyman continuum flux 
inferred from the de-reddened Br7 flux (3.1 x 10~ 14 erg 
s _1 cm -2 ), Q(H + ) = 1.0 x 10 53 photons s -1 , to constrain 
the cluster mass and age. Using instantaneous Starburst99 
models we find a total mass of ~ 7 x 10 6 Mq (with ~ 2600 
O stars) for a Salpeter IMF extending from 1 to 100 Mq, 
and an age of ~ 4 Myr. This age is consistent with the 
lack of photospheric CO and metal absorption lines from 
red supergiants and other cool giants, which would begin 
to contribute significantly to the 2 /im light at an age of 
~ 7 Myr (Leitherer et al. 1999). The cluster's density is 
then about 115 M Q pc~ a for stars of 0.1-100 M Q within 
a half-light radius of 32 pc. This density is 30 times less 
than that of the LMC SSC, R136 (within a radius of 1.7 pc, 
assuming a Salpeter proportion of low-mass stars) (Hunter 
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et al. 1995). Thus the Antennae cluster may be a complex 
of clusters rather than one massive cluster. 



3.1. Nebular Emission 

The cluster spectrum features a variety of nebular lines 
that reveal information about the conditions in and around 
the gas ionized gas by the cluster, which in turn allows us 
to constrain the effective temperature of the ionizing stars. 

We detected H I Pfund series lines from Pf 19 to Pf 38, 
which allow us to infer the extinction across the K win- 
dow. We display the Pfund fluxes relative to that of Br7 
in Figure ||, where filled symbols give fluxes for the blends 
Pf 28+H 2 2-1 S(0) and Pf 29+[Fe in]. They fall well 
above the other points, which follow closely the theoreti- 
cal expectation for intensities relative to Br7 (solid curve) 
with no reddening applied, for a gas with n c = 10 4 cm~ 3 
and T c = 10 4 K ( |Hummcr fe Storey 1987| ). Excluding the 
two known blends, the best-fit foreground screen extinc- 
tion is Ak = 1.1 ± 0.3 mag (dashed curve), assuming the 
extinction law of Landini et al. (1984) and evaluated at 
2.2 /im. We consider this an upper limit on Ak because 
a close look at the spectrum shows that the points above 
the dashed line in Figure | for Pf 22-24 at 2.404, 2.393, 
and 2.383 /im may also be blended or contaminated by 
sky emission, implying a lower Ak and a much better fit 
to the theory. Hence the majority of the extinction to the 
cluster is bypassed by observing it in K band. 

The ratios of [Fe III] 2.146 /zm and 2.243 /im lines to 
[Fe III] 2.218 /im are nebular density diagnostics, consis- 
tcnt with a fairly hig h density, n = 10 35 — 10 4 cm~ 3 . (See 



Gilbert et al. (2000) for more detail.) The He I line ratios 



can be used to infer nebular temperature T c , and are fairly 
insensitive to n . We find the ratio He I 2.1128+2.1137 
/im/He I 2.0589 /im = 0.052 ± 0.003, which f or n c = 10 4 



is consistent with T c = 17, 500 ± 800 K ( [Bcnjamir| 



et al. 199£). However this is much hotter than typical neb- 
ular temperatures, and may indicate some non-nebular 
contribution from hot stars to the line emission. 

The He I 2.0589 /im/Br 7 rati o is an indicator o f the 
T e ff of hot stars in H 11 regions ( Doyon et al. 1992 ), al- 
though it is sensitive to nebular conditions such as the 
relative volumes and ionization fractions of He + and H + , 
geometry, de nsity, dustiness, etc. ( Shields 1995 ). Doherty 



et al. (1995) studied H and He excitation in a sample of 
starburst galaxies and H II regions. For starbursts they 
found evidence for high-T e ff, low-n e (~ 10 2 cm -2 ) ionized 
gas from He 1 2.0589 /im/Br 7 ratios of 0.22 to 0.64. The 
ultra-compact H 11 regions were characterized by higher 
ratios (0.8—0.9) and higher densities, ~ 10 4 cm~ 3 . The 
cluster has a flux ratio of 0.70, a value between the two 



object classes of Doherty et al. (1995). Assuming the line 
emission is purely nebular, this ratio is consistent with a 
high-density (10 4 cm -3 ) model of shields (1993) (also in- 
dicated by the [Fe ill] emission), and implies T c ff ~ 39, 000 
K for the assumed model parameters. 




2.32 2.34 2.36 2.38 2.40 2.42 2.44 
Rest Wavelength (urn) 

Figure 4- Pfund line fluxes relative to Brj flux (1.05 x 
10~ 14 erg s" 1 cm~ 2 ). Solid curve is unextincted theoret- 
ical curve for n c = 10 4 cm -3 , T c — 10 4 K (Hummer & 
Storey 1987). Filled symbols represent lines that are known 
blends, and the dashed curve shows theoretical fluxes with 
the best-fit extinction = 1.1 mag. 



The cluster has properties more like those of a compact 
H 11 region than a diffuse one. It appears to be a young, 
hot, high-density H 11 region, one of the first to form in 
this part of the Antennae interaction region. 

3.2. Molecular Emission 

The spectrum shows evidence for almost pure UV fluo- 
rescence excited by FUV radiation from the O & B stars; 
the strong, vibrationally excited 1—0, 2—1 & 3—2 H 2 emis- 
sion has T vib > 6000 K and T rot ~ 970, 1600, and 1800 
K, respectively, and weak higher-v (6—4, 8—6, 9—7) tran- 
sitions are present as well. The H 2 lines are extended over 
~ 200 pc, about twice the extent of the continuum and 
nebular line emission, so a significant fraction of the FUV 
(912-1108 A) light escapes from the cluster to heat and 
photodissociate the local molecular ISM. 

We compared the photodissoci ation region (PDR) mod- 
els of |Draine &: Bertoldi (1996) with our data by cal- 



culating reduced xt- Models with high densities (njj 
10 5 cm -3 ), moderately warm temperatures (T = 500 to 
1500 K at the cloud surface), and high FUV fields (Go = 
10 3 — 10 5 times the mean interstellar field) can reasonably 
fit the data. Figure |^ shows xt contours for all models 
projected onto the nn— Go plane. The best-fit Draine & 
Bertoldi model is n2023b, which has n H = 10 5 cm~ 3 , T 
= 900 K, and G = 5000. We fit 22 H 2 lines, exclud- 
ing 3—2 S(2) 2.287 /im because it appears to be blended 
with a strong unidentified nebular line at 2.286 /im found 



in higher-resolution spectra of planetary nebulae (Smith 
|st al. 19"81 ). The weak high-v transitions are all under- 



predicted by this model, and appear to come from lower- 
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density gas (nn ^ 10 3 — 10 4 cm 3 ) exposed to a weaker 
FUV field (G £ 10 2 - 10 3 ). 
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Figure 5. Comparison of Hi lint strengths with PDR mod- 
els. Contours of xt f or 2% lines projected onto n^ — Go 
plane peak at na ~ 10 5 cm~ 3 and Go ~ 5000. Model points 
(+) are for T = 300 - 2000 K. White + marks best-fit 
PDR model of Draine & Bertoldi (1996), with T = 900 
K and xl = 9.3. Contours are 50, 25, 20, 15, 12, 10. 



The ortho/para ratio of excited H2 determined from 
the relative column densities inferred from 1—0 S(l) and 
S(0) lines is 1.62±0.07. This is consistent with the ground 
state v=0 H2 being in LTE with ortho/para ratio of 3 if 
the FUV absorption lines populating the non-LTE ex cited 
states are optically thick ( Sternberg fc Ncufcld 199S| ). In- 
deed, the best-fit PDR models have temperatures that are 
comparable with T rot in the lowest excited states. 

If the extent of the H2 emission indicates that the 
mean-free path of a FUV photon is ~ 200 pc, then (nn) 



= 3 cm -3 for a Galactic gas-to-dust ratio, while in the 
PDR(s) nn = 10 4 — 10 6 cm -3 . This implies that the molec- 
ular gas is extremely clumpy, which is consistent with the 
range of densities inferred from the detection of anoma- 
lously strong v = 8—6 H2 emission. 



4. Conclusions 

The highest surface brightness mid-IR peak in the ISO 
map of the Antennae Galaxies is a massive (~ 16xl0 6 
M©), obscured (Ay ~ 9 — 10), young (age ~ 4 Myr) star 
cluster with half-light radius ~ 32 pc, whose strong FUV 
flux excites the surrounding molecular ISM on scales of up 
to 200 pc. The cluster spectrum is dominated by extended 
fluorescently excited H2 emission from clumpy PDRs and 
nebular emission from compact H 11 regions. In contrast, 
the nearby nucleus of NGC 4039 has a strong stellar spec- 
trum dominated by cool stars, where the only emission 



lines are due to shock-excited H2 . These observations con- 
firm the potential of near-infrared spectroscopy for ex- 
ploration and discovery with the new generation of large 
ground-based telescopes. Our ongoing program of NIR- 
SPEC observations promises to reveal a wealth of infor- 
mation on the nature of star formation in star clusters. 
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